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: ABSTRACT 

We investigated the abundance ratio of ^^C/^^C in planetary nebulae by ex- 
amining emission lines arising from C III 2s2p ^-P2,i,o ~^ ^5*0. Spectra were 
•/^ ! retrieved from the International Ultraviolet Explorer archives, and multiple spec- 

^ ! tra of the same object were coadded to achieve improved signal-to-noise. The 

^ , hyperfine structure line at 1909.6 A was detected in NGC 2440. The ^^C/^^C 

Tjj- I ratio was found to be ~4.4±1.2. In all other objects, we provide an upper limit 
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for the flux of the 1910 A line. For 23 of these sources, a lower limit for the 
i2q^i3q j^a^iQ established. The impact on our current understanding of stel- 
lar evolution is discussed. 

The resulting high signal-to-noise C III spectrum helps constrain the atomic 

C/3 
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I physics of the line formation process. Some objects have the measured 1907/1909 

flux ratio outside the low-electron density theoretical limit for ^^C. A mixture of 
^^C with ^^C helps to close the gap somewhat. Nevertheless, some observed 
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^ [ 1907/1909 flux ratios still appear too high to conform to the presently predicted 

limits. It is shown that this limit, as well as the 1910/1909 flux ratio, are predomi- 
nantly influenced by using the standard partitioning among the collision strengths 
for the multiplet ^Sq-^Pj according to the statistical weights. A detailed calcu- 
lation for the flne structure collision strengths between these individual levels 
would be valuable. 

Subject headings: Atomic processes — planetary nebulae: abundances, general, 
individual (NGC 2440) — stars: abundances, AGB and post-AGB 
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1. Introduction 

The relative abundances of ^^C and ^^C in evolved stars present a strong challenge 
to our current understanding of stellar nucleosynthesis. A growing body of observational 
evidence suggests that standard stellar models greatly overestimate the ratio of ^^C to ^'^C 
on the red giant branch (RGB) and later evolution. In a careful study of the literature, 
Charbonnel & do Nascimento (1998) have shown that a full 96% of evolved stars exhibit a 
i2Qy/i3Q ratio below 15, considerably lower than models predict (^^C/^^C between 20 - 30). 
On the other hand, cool bottom processing (CBP), which employs a deep mixing mechanism 
below the standard convective envelope during the red giant phase (Sackmann & Boothroyd 
1999), predicts ^^C/^^C ratios <20 for low-mass post-RGB stars (Boothroyd & Sackmann 
1999). Using millimeter wave observations of ^^CO and ^^CO, Palla et al. (2000) found 
several examples of planetary nebulae (PNs) with small ^^C/^^C ratios, as low as 9 for NGC 
7293. Balser et al. (2002) in a similar study measured ^^CO and ^^CO in molecular clouds 
associated with PNs. They found the range of ^'^C/^^C from 2.2-31 in 9 PNs. The extreme 
low value 2.2±0.03 was found for the PN Ml-16. Developing another approach to obtaining 
these abundance ratios is very valuable to corroborate these low values. 

Clegg (1985) pointed out that the ^^C/^^C ratio could be determined using an extremely 
weak C III ^^C hne at 1909.597 A. This hne is a hyperfine induced transition that arises from 
2s2p ^Pq — > 2s^ ^Sq, which is strictly forbidden in ^^C. By comparing the flux of this hne to 
the nearby bright C III lines 2s2p ^Ps,! ^ 25^ ^Sq at 1906.683 A and 1908.734 A respectively, 
one can find the ^^C/^-'^C ratio. Using the Hubble Space Telescope (HST), Clegg et al. (1997) 
found a ratio of 15±3 in NGC 3918 and 21±11 in SMC N2, as well as a very marginal 7^]^ 
for LMC N122. Palla et al. (2002), also using HST, reported a non-detection in NGC 3242, 
showing ^2c/i3c > 38. 

When studying faint lines in extended sources, the light gathering power of the instru- 
ment is extremely important. The large lUE aperture, an ~20"x~10" 'rounded rectangle' 
(also described as an oval or race-track pattern), allowed a high throughput for extended 
objects like PNs. The extensive archival database and New Spectral Image Processing Sys- 
tem (NEWSIPS) allow us to effectively study the 1909.597 A ^^C hne with lUE. Using these 
data, coadding multiple spectra when available to get a higher signal-to-noise ratio (S/N), we 
attempt to measure or provide upper limits for the flux of this line. In conjunction with our 
measurements of fluxes in the adjacent "main" fines, we derive or set limits for the ^^C/^^C 
ratio in most of the sample of PNs. 

In section 2, we discuss the lUE database and our flux measurements. The analysis tools 
to interpret the observations are introduced in section 3, with the population rate equations. 
In section 4, the flux ratios are analyzed in terms of the ^^C/^^C ratio and the electron 
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density (Ne). Section 5 contains a discussion and conclusions. 

2. lUE Observations 

The data were obtained from the lUE archives. For most sources, only short wavelength 
prime (SWP) camera spectra with the large aperture were used, because the long wavelength 
spectra taken with both the prime camera (LWP) and redundant camera (LWR) are much 
noisier at this wavelength. The spectral resolution of high dispersion SWP, LWP, or LWR 
spectra at these wavelengths is ~0.2 A, sufficient to separate the weak ^^C line from the 
bright hne at 1908.7 A. 

We included all promising PNs observed with high dispersion in the large aperture. 
Multiple observations of the same object were coadded point by point, weighted by the inverse 
of the total noise in the continuum, to improve the S/N, using our own analysis program. 
We note that the inverse noise and exposure time are not linearly related in lUE spectra. 
Figure 1 shows a plot of inverse noise vs. integration time for NGC 3918 SWP spectra. 
Though the data quality, as gauged by inverse noise, generally improves with increased 
time, the quality is not a linear relation of either exposure time or the square root of the 
exposure time. The lUE on-line "Frequently Asked Questions" (http://archive.stsci.edu/cgi- 
bin/faq.cgi?mission=IUE#40) recommends inverse noise weighting, and our data show that 
weighting by the inverse noise generates better S/N than weighting by the time or by the 
square root of the time. 

The ISAP^ routine Line Fit was used to fit a three-component Gaussian to the three 
lines of the C III multiplet in order to determine the line fluxes or to set an upper limit when 
a line was not detected. The number of free parameters were limited by constraining the full 
width half maximum (FWHM) to be that of the 1907 A line. Only one degree of freedom in 
the determination of wavelength was allowed by requiring the fit to conform to the known 
AX^ac between the transitions - that is an offset of 2.051 and 2.914 A for the 1909 and 1910 
lines relative to the 1906.683 A line. Wc use a quadratic continuum fit to the line-free data. 
Most of the lUE spectra that we found useful are SWP spectra. For fitting the continuum 
basehne in these, we avoid the data points at lower wavelength than the 1907 A line where 
there is an echelle interorder "sphce" (at ~1905.4) and the fluctuations in the continuum 
intensity become much larger than those on the red side. 



^ISO Spectral Analysis Package (ISAP) is a joint development by the LWS and SWS Instrument Teams 
and Data Centers. Contributing institutes are CESR, IAS, IPAC, MPE, RAL, & SRON. 
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Fig. 1. — The inverse noise (the reciprocal of the sum of noise points from ~1910. 5-1920 A) 
versus the exposure time for the various SWP high dispersion, large aperture lUE spectra 
of NGC 3918. Units in the ordinate are relative; a larger value indicates less noisy data. 
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We used both Gaussian-fitting and direct integration routines to measure the net hne 
flux above the continuum. Because the fiUed large aperture introduced a spectral impurity 
into the spectrograph, the line profiles have flatter tops and less extended bases (e.g., are 
more "trapezoidal" ) than the Gaussian fits. For this reason and because some of the sources 
have lines that depart in other ways from Gaussian, we prefer to use the direct integration 
routines (the moment option in ISAP) to measure the line flux. We also confirmed our line 
flux results using IRAF^ splot with the e-option for direct integration. 

The root-mean-squarc (rms) fluctuations in the continuum flux were determined from 
a region clear of lines and obvious spikes longward of 1910 A and the FWHM determined 
from the 1906.7 A line. From this, we determine the value 3a, which is our threshold for 
judging a flux measurement a detection or not. This 3a value is also used to set upper 
limits on the flux for the non-detections. The results of our measurements are presented in 
Table 1. Column 1 has the source name; columns 2 and 3 the J2000 positions; columns 4, 6, 
and 8 have the line flux measurements F(1906.7), F(1908.7), and F(1909.6) in erg cm~^ s~^. 
Column 5 has the la uncertainty associated with each of the line fluxes while column 7 is 3 
times this value in order to judge whether or not the hyperfine line has been detected. The 
FWHM from the Gaussian fit of the coadded, unsaturated spectra are in column 9. The 
entries in the remaining columns are quantities derived from these data and will be discussed 
later. 

For many of the sources, there is a second line in Table 1 where more lUE spectra were 
coadded with those used to produce the first-line results (with only unsaturated data) . These 
added spectra are saturated in either the 1907 or 1909 lines. Because of the limited dynamic 
range of lUE, long exposures, particularly with the SWP camera, of bright PNs will readily 
saturate (overexpose) these strong C III lines. However, there is not a saturation issue with 
the 1910 line and adjacent continuum. Using these additional data, we are able to test if 
there are now any statistically significant detections of the 1910 line by comparing F(1909.6) 
with the 3a uncertainty, which is generally smaller than the entry on line 1. Although the 
resulting fiux and FWHM in the main lines are not reliable, we are able to use these lines as 
a fiducial wavelength indicator only to fit the data in the vicinity of the 1910 A line with a 
single Gaussian that has a FWHM assumed to be the same as on line 1. As will be discussed 
(§4.2), the only object with a significant detection of the 1910 hne is NGC 2440. Even with 
the additional data entered on line 2, we found no other significant detections of this line. 



^IRAF is distributed by NOAO, which is operated by AURA, under cooperative agreement with NSF. 
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3. Rate Equations for the Level Populations 

In order to interpret the observational data, we solve the C III statistical equilibrium 
equations for the populations of the 4 lowest-lying energy levels, numbered 1-4. These 
are ^5*0, 2s2p ^Pq, 2s2p ^Pi, and 2s2p ^P2- For this, the effective collision strengths 
(T values) and transition probabilities (A-values) listed in Table 2 are used. Note the 
important distinction for A21 between ^^C and ^^C. This is the sole quantity that is different 
and is the crucial factor in providing different results when solving the rate equations for 
the level populations. The effective collision strengths for transitions between the ^Pj levels 
are available for the electron temperature (Tg) = 5000, 10000, 15000, and 20000 K (Keenan 
et al. 1992). We use these same temperatures in Table 2, adding the recent data for ^5*0- 
^Pj (Mitnik et al. 2003 and Don Griffin, private communication). Then, following standard 
procedure, their value is partitioned among the ^Pj levels according to the statistical weights. 

The results of the solution of the rate equations for ^^C and for ^'^C in terms of the 
pertinent volume emissivity (j value) ratios are given in Figure 2. In this Figure, we plot 
J1907/J1909 for both pure ^^C and pure ^^C in addition to ji9io/ji909 for pure ^^C. These ratios 
are shown vs. the log (iVg cm~^) for two values of Tg, 10000 and 15000 K. This is a typical 
range for Tg in PNs and demonstrates that there is minimal dependence on Tg. 

Besides the graphical representation shown in Figure 2, it is useful to solve the 4-level 
rate equations explicitly in the low-iVg limit in terms of algebraic expressions for the atomic 
data. As will be seen, this provides superior insight when we compare the observed line 
ratios with theoretical predictions. The expression for the familiar A^g-diagnostic line ratio 
may be written, 

J(1909) A^3 ^31 Jysi' ^ ' 

Here we need to obtain A^4/A3. How this was done for the general 3-level atom rate equations 
may be seen from equations 2 and 3 in Rubin (1986). Because we are writing the equations 
for just the low-A'e limit, when ^421 7^ 0, all the terms arising from levels above ground that 
contain the square of the density are dropped (e.g., coUisional routes out are negligible). 
Also terms with very low Aij may be neglected. Then, in the low-Ag hmit when A21 7^ 0: 

^(1^^1001 (2) 

/(1909) ■ Ci3An + CMA43' 

where Cij = Tij exp{—Xij/kTe). The xs are the energy level differences and k is the Boltz- 
mann constant. Because A43 << A41, then 

/(1907) Ci4 
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Table 2. Effective Collision Strengths and Transition Probabilities for 4-Level C III 



Te 



i-j 


5000 


10000 


15000 


20000 


Ajj s 


1-2'' 


0.13 


0.11222 


0.10711 


0.10556 


(for 12 c) 












9.04E-4 (for 13C)<= 


1-3'' 


0.39 


0.33667 


0.32133 


0.31667 


103'''^ 


1-4" 


0.65 


0.56111 


0.53556 


0.52778 


5.14E-3': 


2-3'' 


0.800 


0.950 


1.01 


1.03 


2.39E-7'^ 


2-4^ 


0.590 


0.695 


0.797 


0.855 


1.4E-13'= 


3-4^ 


2.23 


2.75 


3.05 


3.23 


2.39E-6'^ 



''^Don Griffin (private communication) using Mitnik et al. 2003 
partitioning among levels a<;cording to statistical weights 

''Keenan et al. 1992 

"^Brage et al. 1998 

<lDoerfert et al. 1997 

''Jonsson & Frocsc Fischer 1998 

^Fleming et al. 1996 
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When we assume that all exp{—xij / kTe) are approximately the same, then 

7(1907) Ti4 5 



(4) 



7(1909) Ti3 3 ' 

following the ratio of the statistical weights. This is the equation governing the ^^C ratio. 

In the \ow-Ne limit when A21 = 0, the coUisional terms arising from level 2 are retained. 
Thus, in the low-A^e limit when ^421 = 0: 

/(1907) ^ C12 T24 + Cu (T23 + T24 + T12) . . 

7(1909) ■ Cu T23 + Cis (T23 + T24 + T12) ■ ^ ^ 

The last T12 term arises from collisional deexcitation out of level 2 and the symmetry of T21 
and T12. Again, when all exp{—xij/kTe) are nearly the same, then 

/(1907) Tr2 T.i + Tn (T^, + T^, + T12) , . 

7(1909) " T12 T23 + Ti3 (T23 + T24 + T12) ■ ^ ^ 

This is the equation governing the ^^C ratio. We note that the second terms (including the 
parenthesized values) in the numerator and denominator dominate, and if the first terms 
were not present, would give the same result 5/3 as the ^^C ratio. It is interesting that the 
two solutions hinge on just ^421. Yet it is true that the low-A^e limits for both the ^^C and 
^^C C III ion depend on the collision strengths. 

Let us now consider the ratio 7(1910)/7(1909) in the low-A^e hmit with A21 ^ 0. This 
is for the case of pure ^^C. The bottom curve in Figure 2 and the behavior in the low-A^e 
that we are about to work out do not include any contribution from ^^C to the main lines. 
Starting with 

7(1910) ^ N2 A21 V2X 

7(1909) iV3 A31 i/31 ' ^ ^ 

we find, 

/(1910) ^ C7l2 A4, - C,4 ^43 ,^ 

7(1909) ^-^^^^ Ci3 A41 + Ci4 A43 ■ 
The comments and assumptions that lead from equation (1) to (2) apply here, except that 
N2/N^ is apropos now. Again, when all exp{—Xij / kT^) are nearly the same, then 

/(1910) T12 A41 - Ti4 . . 

7(1909) '"Ti3 A41 + T14A43 ■ ^ ' 

The second terms in the numerator and denominator are negligible compared with the first 
terms; by dropping them, one obtains 

7(1910) T12 1 



7(1909) Ti3 3 ' 



(10) 
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Fig. 2. — The ratio of volume emissivities vs. log {Ne cm~^), from top to bottom: jigor/jigog 
for pure ^^C and ^^C respectively, jigio/jigog for pure ^^C. Note that the jigog values differ 
for the ^^C and ^^C curves. The solid and dashed lines are for Tg = 10000 and 15000 K. 
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foUowing the ratio of the statistical weights. 

The equations derived in this section agree extremely well with the low-A^g hmits shown 
in Fig. 2. The small variations with T^. are also in agreement between the analytical expres- 
sions and Fig. 2. 



4. Interpretation of the Observed Line Fluxes 

It is important to realize that the "main" lines 1907 and 1909 are produced by both ^^C 
and -"^^C. Following Brage et al. (1998), we may write an equation for the intrinsic line ratio 
7(1910)//(1909) = hi/ hi as 

hi _ 1^21 ^2^21 

hi ^2.1 {rms + naj/lai 

where r is the ^^C/^^C abundance ratio, n2, n^, and are the fractional populations for 
levels 2, 3, and 4 for ^^C and and 7714 arc the fractional populations for levels 3 and 4 for 
^^C. That is the sum over rii, i = 1-4 is unity, as is the sum over rrij. Values of rii and mj 
are determined from solving the respective rate equations. 

We also write an analogous equation for the ratio of the main hues /(1907)//(1909) = 
hi /hi as 

hi _ i^Ai {rrrii + n4)A4i 
^31 ^zi (rrris + n3)Asi 

For each of the sources in Table 1, there are measured flux ratios for F(1907)/F(1909). In 
addition, an upper hmit to F(1910)/F(1909) may be set by using 3a as the upper limit for 
the non-detected 1910 A lines. (Because lines in this multiplet are so close together, intrinsic 
flux ratios are equivalent to observed flux ratios.) 

Equations (11) and (12) are two relations that depend on r and A^e with only a weak 
dependence on T^. For the solution of these equations, we use — 10000 K. An iterative 
solution converged rapidly with the results for cm"^ and r, the ^^C/^^C abundance ratio, 
entered in the last two columns of Table 1. Just as rapidly, there is an indication when there 
is no possible solution. For these dash is shown in the ^^C/^^C- column with the 

Ne value found assuming pure ^^C. The reason that a range in Ne is shown for those PNs 
with a lower limit for ^^C/^^C is that the higher number for corresponds to the solution 
of equations (11) and (12) for that limiting value of r. This value of N^, is in fact an upper 
limit to the actual N^, which may range down to the smaller number shown for pure ^^C. 
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4.1. The Flux Ratio F(1907)/F(1909) 

For some of the PNs, we found that the measured 1907/1909 A flux ratio was higher 
than the theoretical low-iVg hmit for ^^C 1.511, 1.526, and 1.537 at 10000, 15000, and 
20000 K, respectively. For instance NGC 6818 has F(1907)/F(1909) = 1.73, significantly 
higher than the atomic data predict. Rubin et al. (2001) found a similar problem for other 
A^e-diagnostic fine ratios. In the most recent work, Rubin (2003) found that the observed 
[Ne V] 14.3/24.3 /im flux ratio was out of theoretical bounds on the low-Ag limit for 10 out 
of 20 PNs in the Infrared Space Observatory (ISO) archive, including NGC 6818. 

A mixture of ^^C with ^^C helps to close the gap somewhat. The theoretical low-A^e hmit 
for ^^C is 1.653, 1.658, and 1.660 at 10000, 15000, and 20000 K, respectively. However, 
even if half the carbon were ^^C, observed F(1907)/F(1909) ratios as high as 1.73 would 
not be explained. Indeed any ratio exceeding the 5/3 value would not comport with present 
predictions. The value of 5/3 is the result of the standard practice of partitioning the ^Sq-^Pj 
multiplet value among the flne-structure levels according to the statistical weights (1, 3, 5). 
Consequently, this 5/3 limit does not depend on whose set of collision strengths is employed 
as long as they are partitioned in this way. An observational program to measure very 
accurate F(1907)/F(1909) ratios would be most useful. On the theoretical side, a detailed 
calculation for the effective collision strengths between these individual LS levels would be 
extremely valuable. 



4.2. The Flux Ratio F(1910)/F(1909)- Detection of the line in NGC 2440 

The lUE data and analysis presented in Table 1 show only one source, NGC 2440, 
with a signiflcant detection of the ^^C line (Figure 3). This is the coadded spectrum of two 
archival spectra SWP 7263 and LWR 10741. The apertures for both spectra were centered 
at the same position. The position angles of the aperture major axis were 199.19 and 20.32° 
respectively (only ~1° mismatch) thus providing nearly the same sky coverage. Exposure 
times were 3600 s on 1979 November 29 and 16200 s on 1981 May 29 respectively and were 
observed under programs NCBJH and NPDJH, PI J. Patrick Harrington. 

In the coadded spectrum the observed flux of the 1909.6 A line is (5.94±1.47) xl0~^^ 
erg cm~^ s~^, which is a 4<t detection. The simultaneous solution of equations (11) and (12), 
assuming = 10000 K, gives A^e = 2460 cm"^ and ^^C/^^C = 4.4. With an assumed Tg of 
15000 K, the solution becomes = 3130 cm~^ and ^^C/^^C = 4.2. We prefer the former 
solution because according to figure 2 (Tg vs. LP.) in Bernard Salas et al. (2002), it appears 
that the Te value for the C++ zone in NGC 2440 is closer to 10000 than 15000 K for the case 
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1906 1908 1910 1912 

Wavelength (A) 



Fig. 3. — Spectrum in the vicinity of the C III multiplet for NGC 2440 resulting from the 
coaddition of lUE high dispersion, large aperture spectra SWP 7263 and LWR 10741. The 
continuum is very weak compared with the 1907 and 1909 A lines and has been fit with a 
quadratic baseline and removed. The inset shows a magnified (in intensity) view of the ^"^C 
1910 A line. The 3-component Gaussian fit (see text) and sum are shown, as well as the line 
connecting the data points (diamonds). 
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0+ — > C''""'" with ionization potential (LP.) 24.4 eV. Although not sensitive to the assumed 
Tg, the general trend we find here is that the higher the Te, the higher the derived iVg and 
the lower the derived ^^C/^^C. We estimate the uncertainty in the value of ^^C/^^C as ±1.2 
from separate measurements and analysis of the two individual lUE spectra. 

In their study of NGC 2440, Bernard Salas et al. used these same two lUE spectra. Ac- 
cording to their table 3, F(1907) = 2.43x10"" and F(1909) = 1.77x10"" erg cm"^ s-\ This 
leads to F(1907)/F(1909) = 1.373, while they have a slightly different ratio of 1.35 in their 
table 5 where they also list the derived = 4360 cm""^ assuming = 15000 K. We compare 
these to our entries in Table 1: F(1907) = 2.24x10"" and F(1909) = 1.57x10"" ergcm'^s'^ 
This leads to F(1907)/F(1909) = 1.427, which will result in a lower derived A^e irrespective 
of accounting for the effect of the presence of ^^C. It is not clear why their hue fluxes are 
higher than ours. It may be due to several factors. For instance, they may have used the pre- 
NEWSIPS cahbration. They likely used a different technique than ours to form a weighted 
average of the two lUE spectra. Their fluxes may be taken directly from the Gaussian-fit 
areas whereas we have used Gaussian fits for guidance but use the option in ISAP [moment 
fit] and IRAF [splot e-option] to measure the areas (except for measuring the 1910 line with 
the Gaussian). Generally we find that departures from a Gaussian profile result in some- 
what smaller line fluxes. Additionally, the fltting of lines is not totally objective, e.g, in the 
manner of accounting for the baseline, etc. 

Although the line was not detected in any of the other PNs, we were able to establish 
lower limits for the ^^C/^^C ratio for many. Even when coadding many spectra, the S/N 
was not high enough to ascertain a useful value of ^^C/^^C for most objects. HST studies 
similar to those of Palla et al. (2002) and Clegg et al. (1997) could be used to refine the 
limits presented here. 

At a late stage in our work on this project, we were surprised to discover that a similar 
project had been undertaken by Miskey et al. (2000). They report a detection in NGC 2440, 
but appear to derive a much higher ^^C/^^C ratio of 39. They also report a detection in 
NGC 6302 with a ratio of 23. Unfortunately, since only an AAS abstract has been pubhshed, 
their data and methods are not available for scrutiny. We cannot report a detection of the 
1910 A line in NGC 6302. Two short wavelength spectra (SWP 10391 and SWP 8971) and 
one long wavelength spectrum (LWR 7722) are available for this nebula. Even the 1907, 
1909 main lines are indistinguishable from the noise in LWR 7722; coadding SWP 10391 and 
SWP 8971 produces a non-detection of the 1910 line. SWP 10391 shows a small peak at the 
correct location, but at the same level as the noise; thus a detection cannot be conclusively 
claimed. 



While ten SWP and eight LWP/LWR spectra were available for NGC 3918, many of 
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them showed saturation in the strong C III hnes, making them unusable for determining 
ratios. Similarly, all four available SWP spectra for NGC 3242 are saturated, and the values 
from Tabic 1 come from only LWR 12705. Our limits are consistent with the findings of 
Clegg ct al. (1997) and Palla et al. (2002). We note that Palla et al. used a la upper limit 
for F(1910), while we use 3a. Using their la value and main line fluxes for NGC 3242 with 
our improved atomic data here, we find that their ^^C/^^C > 38 becomes > 44.4. However, 
when we use a 3a observed upper limit from their paper, then more conservatively ^^C/^^C > 
14, a value very similar to our limit in Table 1. 



5. Discussion and Conclusions 

NGC 2440 has an unusually massive central star (CS), 0.72 M0, according to Wolff ct al. 
(2000), and 0.932 M©, according to Corny et al. (1997). The progenitor star is believed to be 
more massive than 4 Mq, e.g., figure 3 in Balser et al. (2002) has ~4.2-4.3 Mq (see also the 
discussion in their §2). In their study of NGC 2440, Bernard Salas et al. (2002) concluded, 
based on a measured (C+N+0) /H ratio higher than that for NGC 7027, that the CS must 
have been a more massive progenitor than was the star in NGC 7027, for which Bernard 
Salas et al. (2001) found 3-4 Mq. Above 4 Mq, asymptotic giant branch (AGB) stars develop 
deep convective envelopes with very high base temperatures, leading to the operation of the 
CN-cyclc; hence the name, hot bottom burning. There, the CN-cycle operates at the base 
of the convective envelope in a hydrogen-burning shell, where ^^C is converted into ^^C and 
then into ^^N. This process would tend to drive ^^C/^^C toward the CN-cycle equilibrium 
ratio of ~3.3 (Smith & Lambert 1990; Frost et al. 1998). When ratios close to the CN-cycle 
equilibrium value are measured in a star's atmosphere, it is indicative of the presence of 
nearly pure CN-cycle material. Thus, the low ^^C/^^C ratio of 4.4±1.2 that we find for 
NGC 2440 does support this picture. 

If the progenitor were less massive than 4 Mq, then the scenario changes. Both the cool 
bottom processing (GBP) and standard models predict ^^C/^^C ~ 20-30 for stars between 
3-4 Mq (Boothroyd & Sackmann 1999; Renzini & Voh 1981). Then our derived ^^C/^^C 
for NGC 2440 would not appear to support either the standard or GBP theories. While 
CBP does produce ratios as low as 4.4, it is only expected to occur in progenitor stars of 
<2 Mq (Wasserburg et al. 1995). Two (unlikely) possibilities might allow agreement. First, 
were the mass estimates for the CS progenitor too large by a factor of ~2, our ^'^C/^'^'C = 
4.4 for NGC 2440 would support the CBP theory. Second, perhaps CBP occurs in stars of 
3-4 Mq. This is in fact contrary to the statement by Wasserburg et al. (1995) that no CBP 
is expected above ~2.3 Mq. However a reconsideration of model parameters might be useful 
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to determine whether this could occur. 

The small ^^C/i^C we find for NGC 2440 is in the range derived by Balser et al. (2002) 
in their study of molecular clouds associated with PNs. As mentioned in the Introduction, 
from observations of ^^CO and ^^CO millimeter lines, they found values of ^^C/^^C from 
2.2-31 in 9 PNs. The extreme low value 2.2±0.03 was for the bipolar FN Ml-16, which is 
not on our list. We note that our small ratio for NGC 2440 is consistent with measurements 
of J-type carbon stars. Abia & Isern (1997) found <^^C/^^C> = 6±3 while Ohnaka & Tsuji 
(1999) in their analysis of 26 J-type C stars found an even lower average of 4.7±2.8 with 
several ratios as low as ~l-2, well below the equilibrium value for the CN-cycle. 

The question remains as to why NGC 2440 is the only nebula among those we measured 
which has a detectable 1909.6 A ^^C hyperfine structure line, as well as why it has such a 
low ^^C/^'^C. Using lUE archival data, we also set upper limits on F(1909.6) in 40 other 
PNs available and lower limits for ^^C/^^C in 23 of these. lUE had limited dynamic range 
that precludes setting more stringent limits for most objects, though coadding spectra allows 
better S/N. NGC 2440 is among the brightest PNs in terms of either of the main line fiuxes; 
only NGC 2867, NGC 3242, NGC 3918, and NGC 7662 are higher (see Table 1). By using 
the measured 3a as an upper hmit to their F(1909.6), we derive lower limits on ^^C/^^C for 
these four nebulae of respectively 5.9, 14.5, 9.9, and 6.5 - already higher than what we find 
for 2440. NGC 2440 is carbon rich with C/0 of 1.9 (Bernard Salas et al. 2002) and a high 
excitation Peimbert type I PN (Hyung & Aller 1998), that is, it is a N- and He-rich PN. 
The central star is among the hottest (~200,000 K) and most massive known (see earlier 
discussion). In fact, in their entire sample of ~80 PNs, Gorny et al. (1997) derived a CS 
mass of 0.932 M©, which was the second largest of all. This nebula also has reaching 
as high as 15000 - 20000 K as measured in several high ionization species (Bernard Salas 
et al. 2002). NGC 2440 has a relatively low of 2460 cm'^ determined here (see Table 1). 
Both lower A^^ and lower ^^c/i^c lead to a higher F(1910)/F(1909). It may be that Type I 
PNs have lower ^^C/^'^C ratios. Furthermore, NGC 2440 is a bipolar PN. In their study, 
Balser et al. (2002) found those PNs that have the lowest values of ^^C/^^C are classified 
as bipolar. Further theoretical and observational work is needed to determine the cause of 
the low ratios. An HST spectroscopic study of the C III multiplet in NGC 2440 would be 
desirable. 

Several of the PNs have the measured F(1907)/F(1909) outside the low-Ag theoretical 
limit for ^^C. Uncertainties in the measured fiuxes may account for some of these differences. 
A mixture of ^^C with ^^C helps narrow the gap somewhat. Nevertheless, some of observed 
F(1907) /F(1909) ratios still appear too high to conform to the presently predicted hmits. We 
have shown that both the F(1907)/F(1909) and the F(1910)/F(1909) ratios in the low-A^e 
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limit are predominantly influenced by using the standard partitioning among the collision 

strengths for the multiplet ^Sq-^Pj according to the statistical weights. A detailed calcula- 
tion for the effective collision strengths between these individual levels would be valuable. 
The determination of reliable Ng values by this method as well as the ^^C/^^C abundance 
ratio depends on it. As far as we know, it appears that there have been no previous attempts 
to account for the presence and effects of ^^C when using the well-known C III iVg-diagnostic 
F(1907)/F(1909). 
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